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Abstract 



Plasmodium of Physarum polycephalum is a single cell visible by unaided eye. Dur- 
ing its foraging behavior the cell spans spatially distributed sources of nutrients 
with a protoplasmic network. Geometrical structure of the protoplasmic networks 
allows the Plasmodium to optimize transport of nutrients between remote parts of 
its body. Assuming major Mexican cities are sources of nutrients how much struc- 
ture of Physarum protoplasmic network correspond to structure of Mexican Federal 
highway network? To find an answer undertook a series of laboratory experiments 
with living Physarum polycephalum. We represent geographical locations of major 
cities (nineteen locations) by oat flakes, place a piece of Plasmodium in Mexico city 
area, record the Plasmodium's foraging behavior and extract topology of nutrient 
transport networks. Results of our experiments show that the protoplasmic network 
formed by Physarum is isomorphic, subject to limitations imposed, to a network of 
principle highways. Ideas and results of the paper may contribute towards future 
developments in bio-inspired road planning. 
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1 Introduction 



Plasmodium of Physarum polycephalurr\ ^ \ is a single cell with many diploid 
nuclei. The Plasmodium feeds on microbial creatures and microscopic food 
particles. When colonizing its habitat the Plasmodium develops an optimal 
network of protoplasmic tubes or veins. The protoplasmic network is optimal 
in a sense that it maximizes number of food sources occupied, minimizes time 
for nutrient transport between distant parts of the Plasmodium's body, and 
maximizes area of substrate covered by sensorial activity. 

A Plasmodium is a single cell yet can be considered as large-scale collective 
of simple entities with distribute and massive-parallel sensing, computation 
and actuation. A sensing is parallel because the Plasmodium can detect and 
determine position of many source of chemo-attractants, including nutrients, 
and also make decentralized sensing of environmental conditions, including 
humidity, temperature and illumination. An actuation is parallel because the 
Plasmodium can propagate in several directions in parallel, the Plasmodium 
can occupy and colonize many food sources at the same. With regards to paral- 
lel computation, the Plasmodium is a wave-based massively-parallel reaction- 
diffusion chemical computer pTTlfTO j. A computation in the Plasmodium is 
implemented by interacting bio-chemical and excitation waves [19], redistribu- 
tion of electrical charges on Plasmodium's membrane [Ij and spatio-temporal 
dynamics of mechanical waves [19j. Experimental proofs of P. polycephalum 
computational abilities include approximation of shortest path [22j and hier- 
archies of planar proximity graphs [6] , computation of plane tessellations |31j , 
implementation of primitive memory [28] , execution of basic logical computing 
schemes [36j, control of robot navigation [37j, and natural implementation of 
spatial logic and process algebra [30]. See overview of Physarum-based com- 
puters in [To] . 

Approximation of shortest or even computation of a transportation network 
have been already a hot application for unconventional computing scientists. 
Nature-inspired computing paradigms and experimental implementations were 
successfully applied to calculation of a minimal-distance path between two 
given points in a space or a road network. Thus, a shortest-path problem is 
solved in experimental reaction-diffusion chemical systems [2], gas-discharge 
analog systems [27j, spatially extended crystallization systems [8j, formation 
of fungi mycelian networks [15], and using computer and mathematical mod- 
els of collective insects [12] and Physarum polycephalum [3l]. Previously [5] 
we have evaluated a road-modeling potential of P. polycephalum^ however no 
conclusive results were presented back in 2007. A step forward biological- 
approximation, or evaluation, of man-made road networks was done in our 
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previous paper on approximation of United Kingdom motorways by Plasmod- 
ium of P. polycephalum [9]. Results obtained in [9j were refreshing yet convinc- 
ing. In general, we shown that transportation links constructed by Physarum 
match man-made motorways, there are however subtle differences between 
slime mould and motorways, see details in [lOj. The experimental laboratory 
approach to slime-mould-based approximation of road networks was tested on 
just one country (United Kingdom) so far. To make this bio-inspired approach 
truly universal we must study how the Plasmodium will behave on different 
setups of cities and different shapes of countries. In present paper we have 
chosen Mexico as a test country for Plasmodium approximation because it is 
the most populated Spanish-speaking country with diverse landscapes, sudden 
north-south changes in land size, highest level of public transport in Latin 
America, and high-concentration pollutants in central part of the country. 
These features are not unique but enhanced and, what it is most important, 
they provide very good test-bed for studies of nature-inspired approaches to 
road planning. 

Mexico ranks 7th in highways length with 6,335 kilometers (3,935 miles). Fol- 
lowed by (1) United States with 91,541 km (56,859 mi), (2) China with 24,474 
km (15,202 mi), (3) Germany with 11,400 km (7,081 mi), (4) France with 
10,300 km (6,398 mi), (5) Spain with 9,063 km (5,629 mi), and Italy with 
8,957 km (5,563 mi) p] Political division constitute 31 states and one Federal 
District (DF or Mexico city) that centralize the main activities in the whole 
republic, for this reason main highways at least connect one time Mexico city 
from any place. With a total population of 103,263,388 (2005 estimation) on 
a territorial extension of 1,964,375 km^ (758,450 mi^)[^ Mexico has a relevant 
and strategical geographical position because it is a natural land transit be- 
tween North America (mainly United States and Canada) with Mexico self. 
Central and South America, i.e., all Latin America countries. This contributes 
towards essential components of economical structure for the United States 
in transporting prime materials, natural resources, imports and exports, and 
labor migration. For this reason years ago extensive highways were built from 
Central America across Mexico to the United States, and from coast to coast. 

The paper is structured as follows. In Sect. [2] we describe experimental setup 
and software analysis. Experimental laboratory results on Physarum approxi- 
mation of transport links are presented in Sect. |3) Analysis of Physarum and 
man-made highways in terms of planar proximity graphs is provided in Sect.|4} 
Directions of future studies are outlined in Sect. [5l 



^ 'http : //www . publicpurpose . com/hwy-worldmotorway . htm' 

^ National Institute of Statistical and Geographic (INEGI Spanish abbreviation) 
[http : //inegi . org . mx/^ 
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(a) 




Fig. 1. Experimental basics, (a) Contour map of Mexico with 19 sources of nutrients 
indicated, (b)-(c) Snapshots of two typical setups: urban areas are represented by 
oat flakes, Plasmodium is inoculated in Mexico city, the Plasmodium spans oat flakes 
by protoplasmic transport network. Note that conflgurations of protoplasmic tubes 
are different in these two experiments. 

2 Methods 



Plasmodium of P. polycephalum is cultivated in plastic container, on paper 
kitchen towels sprinkled with still drinking water and fed with oat flakeQ 
For experiments we use 120 x 120 mm polystyrene square and round Petri 
dishes. We use 2% agar gel (Select agar, Sigma Aldrich). Agar plates are cut 
in a shape of Mexico. 

We selected 19 most populous urban areas, based on a selection of relevant 
^ Asda's Smart Price Porridge Oats 
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Fig. 2. General data regions U (Fig.[l|i), including main town and state, population, 
access to international air and sea ports, and economic potential (INEGI source 
2010). 



features, including economic impact, sea and air ports, and tourist attractions. 
Details of the areas are provided in Fig. [2} 

Further we refer to the urban regions as U. These region of U are projected 
onto gel and oat flakes, of size approximately matching size and shape of the 
regions, are placed in the positions of the regions (Fig. ^p). At the beginning 
of each experiment a piece of Plasmodium, usually already attached to an 
oat flake, is placed in the region corresponding to Edo. Mexico, DF, Puebla 
(region 14 in Fig. [T^). 

The Petri dishes with Plasmodium are kept in darkness, at temperature 22- 
25 C^, except for observation and image recording. We undertook 26 experi- 
ments. Periodically the dishes are scanned using Epson Perfection 4490 scan- 
ner. Scanned images of dishes are enhanced for higher visibility, saturation 
increased to 55, and contrast to 40. 

To ease readability of experimental images we provide a binary version of 
each image. The binarization is done as follows. Each pixel of a color image is 
assigned black color if red R and green G components of its RGB color exceed 
some specifled thresholds, R > 9r^ G > Oq and blue component B does not 
exceed some threshold value B < 9b] otherwise, the pixel is assigned white 
color (exact values of the thresholds are indicated in the flgure captions as 
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3 Transport links via foraging 



Snapshots of Plasmodium foraging patterns during typical experiment are 
shown in Fig. [3| Initially a piece of Plasmodium is placed onto an oat flake 
representing Mexico City. The Plasmodium grows and propagates in south- 
easterly direction occupying Oaxaca-Huatulco region and then traveling to- 
wards Tuxtla Gutierrez region (Fig. |3^f). The Plasmodium reaches Tuxtla 
Gutierrez region usually in 12 hours after inoculation (real time distance is 
11 hours). At the same time the Plasmodium colonizes Xalapa- Veracruz and 
Chilpancingo-Acapulco regions and moves north-north- west. In a normal situ- 
ation (vigorous Plasmodium, fresh agar gel, air not contaminated by bacteria), 
south-east and north-north- west developments occurs simultaneously. In 12- 
14 hours after inoculation the Plasmodium can reach as far as Monterrey and 
Nuevo Laredo regions (Fig.|3j3g). The Plasmodium colonizes almost the whole 
Mexico in 38 hours, spanning major urban region from Tijuana in the north- 
west to Tuxtla Gutierrez in the south-east (Fig.jsj^h). The colonization is truly 
completed in 52 hours, when an oat flake representing Merida-Cancun is cov- 
ered by Plasmodium mass (Fig. |3]li). Despite spanning all cities marked by 
oat flakes the Plasmodium continues exploring the territory outlined by agar 
gel and propagates along Baja California Peninsula (Fig. [3^j). There are no 
cities represented by food sources, therefore the Plasmodium retreats from the 
Peninsula. 



In some situations the Plasmodium does not span all cities, represented by oat 
flakes. An example is shown in (Fig. [4]). In few hours after being inoculated in 
Mexico city the Plasmodium propagates in all directions. In 12 hours the Plas- 
modium reaches Ciudad Victoria region on the north, and Xalapa- Veracruz on 
the south (Fig. [4^e). On 24th hour of development the Plasmodium expands 
its occupation till Nuevo Laredo region on the north and Tuxtla Gutierrez 
on the south (Fig. |4|3f). The Plasmodium continues then to Chihuahua and 
Merida-Cancun regions (Fig. ^g) and abandons the oat flake, representing 
Nuevo Laredo region. Even when the whole territory of Mexico becomes col- 
onized by Plasmodium the Nuevo Laredo regions remains free (Fig. lil: 



The Plasmodium may not stop its foraging activity even when all sources of 
nutrients are occupied and the whole agar plate is explored. As illustrated in 
Fig. [5] a vigorous Plasmodium can spread over surrounding Petri dishes, trying 
to settle on a bare plastic. Plasmodium parts residing on a non-agar substrate 
suflFer from the lack of humidity and therefore cease to sustain. Abandoned 
protoplasmic tubes are clearly visible in southern half of Petri dish in Fig.[5]lh. 



As illustrated in Figs. |3pl and [5] Plasmodium rarely develops exactly the 



same foraging pattern twice. Even in any given experiment the Plasmodium 
sometime may change topology of its protoplasmic networks, abandon and 
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(a) t =12 h 



(b) t =24 h 




(c) t =38 h (d) t =52 h 




(e) t =61 h 



Fig. 3. Typical Plasmodium development: (a)-(e) scanned image of experimental 
Petri dish. Time elapsed from inoculation is shown in the sub-figure captions. 
(f)-(j) binarized images, O = (100, 100, 100). 
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(f) t =12 h (g) t =24 h 




(h) t =38 h (i) t =52 h 




(j) t =61 h 
Fig. 3. Continued. 

re-colonize sources of food. Thus we could not ever consider a stationary con- 
figuration of a protoplasmic network but a probabilistic graph, representing all 
possible configuration of protoplasmic networks occurring in experiments for a 
given configuration of sources of nutrients. We define a probabilistic Physarum 
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(a) t =12 h 



(b) t =24 h 




(c) t =35 h (d) t =47 h 



Fig. 4. Plasmodium does not always span all cities (sources of food): (a)-(d) scanned 
image of experimental Petri dish. Time elapsed from inoculation is shown in the 
sub- figure captions, (e)-(h) binarized images, O = (100, 100, 100). 

graph in the following manner. 

Physarum graph is a tuple P = (U, E, w)^ where U is a set of nineteen urban 
areas, E is a set edges, and i(; : E ^ [0, 1] is a probability-weights of edges from 
E. For every two regions a and b from U there is an edge connected a and b if 
a Plasmodium's protoplasmic link is recorded at least in one of k experiments, 
and the edge (ab) has a probability calculated as a ratio of experiments where 
protoplasmic link (ab) occurred to the total number of experiments k. We do 
not take into account exact configuration of the protoplasmic tubes but merely 
their existence. We also consider threshold 9 G [0, 1] Physarum graphs P(^), 
defined as follows: for a, 6 G U (ab) G E if w{ab) > 9. 

Phys arum- gidiphs extracted from 26 laboratory experiments are shown in 
Fig. |6| The graph becomes planar when we remove edges with weights be- 
low 0.19 (Fig. The graphs is acyclic, or a tree, when only edges with 
probability- weights exceeding 0.58 are taken into consideration (Fig.[6]l). Thus 
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(e) t =12 h 



(f) t =24 h 




(g) t =35 h (h) t =47 h 

Fig. 4. Continued. 



edges of the spanning tree are represented by protoplasmic tubes in over half 
of experimental trials. 

Finding 1 Spanning tree is a stable core structure of Physarum foraging and 
nutrient transport network built on a configuration of urban regions U. 

If we increase threshold value to 0.61 the Physarum transport graph becomes 
disconnected (Fig. [7^). The nodes corresponding to Tijuana and Mazatlan 
regions becomes isolated. There is no a direct highway between such cities, 
mainly with main big highways as M15 or M2 and this isolation of course, 
represents geographical limitations. 

The graph is split into two connected components when 9 = 0.65 (Fig. ^]p). 
Superior component reflects directly the highway dominating north-west of 
Mexico that is M2 crossing from Tijuana to Ciudad Juarez. The component 
also represents the fact that Indeed M2 does not reach east coast, that re- 
flect precisely this component. Second component display a strong connection 
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(a) t =24 h 



(b) t =36 h 




(c) t =48 h (d) t =70 h 



Fig. 5. Plasmodium spreads beyond 'dedicated' experimental domain: 
(a)-(d) scanned image of experimental Petri dish. Time elapsed from inoculation 
is shown in the sub-figure captions, (e)-(h) binarized images, O = (100, 100, 100). 

between Mexico city and south-east of Mexico. Such relation involves a high- 
through connection region of north-east with Mexico city, and Mexico city 
with south-east, they are connected specifically by motorway M85 from Nuevo 
Laredo to Mexico city and M180 that run from Xalapa- Veracruz to Cancun. 

A specific relation of connections by standard highways (not rural ways) is 
given in the tablejsjgiven from the last update from "Guia Roji" publications p] 

Finding 2 Transport links to Tijuana and Mazatldn urban regions are unsta- 
ble. 

As we saw some implication in Finding 1 (first component), and table in Fig.js] 
connectivity. We will need at least two or three highways to reach Mazatlan 

^ Source: Guia Roji "Por las Carreteras de Mexico 2011", 17a. Edition. Web site: 
: / / www . guiaro j i . com . mx 



11 



(e) t =24 h 



(f) t =36 h 




(g) t =48 h (h) t =70 h 

Fig. 5. Continued. 



from Tijuana. Here there is not a direct motorway as could be M15. Although 
there is a number of connections to travel between both cities. 

On the way, how well Physarum graphs approximate Mexican highway net- 
work. A sketch of the highway network and the graph H derived are shown 
in Fig. [9} We construct the highway graph H as follows. Let U be a set of 
nineteen urban regions, for any two regions a and b from U, the nodes a and 
b are connected by an edge (ab) if there is a highway starting in vicinity of 
a and passing in vicinity of b and not passing in vicinity of any other urban 
area c G U. 

Finding 3 Graph P(0.58) is a sub-graph o/H apart of three edges. 

Intersection P H H of Physarum and highway graphs is shown in Fig. [9]3. The 
three edges of P(0.58) missing in H are transport links connecting Chihuahua 
and Nuevo Laredo regions, Xalapa- Veracruz and Oaxaca-Huatulco regions, 
and Tuxtla Gutierrez and Merida-Cancun regions. Such difference could pos- 
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(g) 

Fig. 6. Configurations of thresfiold Physarum-graph. F{6) for (a) = 0, (h) 6 = 0.19, 
(c) e = 0.38, (d) 9 = 0.33, (e) = 0.5, (f) 9 = 0.54, (g) 8 = 0.58. Thickness of each 
edge is proportional to the edge's weight. 
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(a) 



(b) 



Fig. 7. Configurations of threshold Physarum-giSiph P{0) for (a) 6 = 0.61 and 
(b) 6 = 0.65. Thickness of each edge is proportional to the edge's weight. 
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Fig. 8. Main highways connections between cities of U. 

sibly be explained by the fact that in laboratory experiments we did not 
represent geographical landscape in agar gel structure but just a uniform gel 
plate was cut in a shape of Mexico. 

Indeed, we have considered natural restrictions. Missing connections between 
Chihuahua and Nuevo Laredo are limited by Chihuahua desert. While con- 
nection between Xalapa- Veracruz and Oaxaca-Huatulco mainly is limited by 
geographical features such as Juarez mountain range. So missing connection 
between Tuxtla Gutierrez and Merida-Cancun implies a large border over Golf 
of Mexico upping to motorway Ml 80 and crossing also on Chiapas mountain 
range. 



4 Physarum, proximity graphs and highways 

Physarum constructs planar proximity graphs by its protoplasmic network [6] . 
A planar graph consists of nodes which are points of Euclidean plane and edges 
which are straight segments connecting the points. A planar proximity graph 
is a planar graph where two points are connected by an edge if they are close 
in some sense. Usually a pair of points is assigned certain neighborhood, and 
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(c) 

Fig. 9. Highway network is highlighted in (a) and graph H of man-made motorway 
network is shown in (b). Intersection P P| H of Physarum P and highway H graphs 
is shown in (c). 

points of the pair are connected by an edge if their neighborhood is empty. 
Relative neighborhood graph [H], Gabriel graph [24], ;5-skeletons [TT] and 
spanning tree are most known examples of proximity graphs. 

Relative neighborhood graph [35] RNG (Fig. flOk) and Gabriel graph [T31I24] 



GG (Fig. 10 3) are computed over nodes corresponding to urban areas U. 
Points a and b are connected by an edge in RNG if no other point c is 
closer to a and b than dist{a^ b) [35j. Points a and b are connected by edge in 
GG if a disc with diameter dist{a^ b) centered in middle of the segment ab is 
empty [T31I24] . The graphs are related as RNG C GG pSllMTi] . 

Proximity graphs found their applications in many fields of science and engi- 
neering. The application domains related to topics of the paper include rout- 
ing in ad hoc wireless networks [T81I32I29II25II38 ]. road network analysis 



study percolation [llj and analysis of magnetic field [33j. Structure of prox- 
imity graphs represents so wide range of natural and engineering system that 
it would be productive to compare Physarum graph P and highway graph H 
to relative neighborhood graph RNG, /5-skeleton BS(1.5) (/3-skeleton B was 
calculated for ^ = 1.5 to make an 'intermediate' graph between RNG and 
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(a) RNG 



(b) BS(1.5) 




(c) GG (d) ST 

Fig. 10. Proximity graphs constructed on regions U: (a) relative neighborhood graph 
RNG, (b) /3-skeleton with control parameter 1.5, (c) Gabriel graph GG, (d) Min- 
imum spanning tree ST 



GG), Gabriel graph GG and minimum spanning tree MST (Fig. 11). From 
now on we mean P(0.19) when writing P. 

Finding 4 RNG C P 



This is a direct outcome of calculating intersections of the graphs (Fig. 11). 
Relative neighborhood graph is considered to be optimal in terms of total 
edge length and travel distance, and is a closest approximation of road net- 
works [391140] . The fact this graph is a sub-graph of Physarum graph indi- 
cates deep intrinsic analogies between formation of man-made and Plasmod- 
ium transport networks. 



In all cases (Fig. lli-d) the large pathway between Mexico city and USA 
is preserved in three following ways. Mexico city via Ciudad Juarez by mo- 
torway M45, Mexico city via Nogales by motorway M15, and Mexico city 
via Nuevo Laredo by motorway M85. They are very important since Mexico 
colonization because they represent ancient ways to transport materials and 
nature sources to Spanish monarchs and later British colonies (now USA). 
Of course, they were besides guided initially from Xalapa- Veracruz to Mexico 
city (from Aztec's age where Spanish tried to reach Tenochtitlan) , and con- 
sequently defining this way from old times, that actually represent motorway 
M150. 
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(a) PflR-NG 



(b)PnBS(1.5) 




(c)PnGG (d)PnMST 




(e) H n RNG (f) H fl BS(1.5) 




(g)HnGG (h)HnMST 



Fig. 11. Intersection of Physarum graph P (a)-(d) and highway graph H (e)-(h) 
with proximity relative neighborhood graph RNG, /3-skeleton RNG(1.5), Gabriel 
graph GG and minimum spanning tree MST. 
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In this context it is not surprising that "Ciudad Juarez" play a strategic posi- 
tion in reaching USA directly from Mexico center and beyond. This is because 
motorway M45 is connected with Ml 90 and runs across all Mexican Republic 
south-to-north from Central America. It is also important to highlight that 
Ciudad Juarez undergone an explosive growth in last few years, partly to its 
dramatic life at the intersection of main drug transport arteries. 

South-east intersections represent connectivity with highways M180 and M190. 
Both are the main large highways connecting principal cities in their regions 
with Mexico city. 



Finding 5 MST C P 



The result is expected because spanning tree is a sub-graph of a relative 
neighborhood graph [35j, however we intentionally highlighted it to show 
that — in terms of minimal-length travel — graph, built by Plasmodium of 
Physarum polycephalum offers optimal solution for transportation of nutri- 
ents. Experiment- wise we found that the spanning tree MST is a sub-graph 



of Physarum graphs P(^) for 9 < 0.58, compare Fig. p^and 111 



As we demonstrated previously (Finding 1) by cutting edges of P(^) with 
9 < 0.58 we transform cyclic graph to a spanning tree (Fig. [6^). By compar- 
ing Fig. [6^ and Fig. [TOl d we find that if Physarum built a transport link con- 
necting Mazatlan and Guadalajara regions instead of connecting Chihuahua 
and Nuevo Laredo regions then the Physarum would approximate an ideal 
minimum spanning tree. Actual travel distance makes Physarum spanning 
tree P(0.58) just a bit longer than ideal spanning tree MST, constructed by 



conventional algorithm Fig. 101. 



In contrast to Physarum graphs highway graph H poorly matches proximity 
graphs (Fig. [TT^-h). 



Finding 6 Let G e {RNG, BS(1.5), GG, MST} then GRH Z5 a discon- 
nected graph. 



Merida-Cancun region stays isolated in all four intersections (Fig. [TTp-h) of 
the highway graph H with proximity graphs. A connected cluster of regions 
Nuevo Laredo, Monterrey, Ciudad Victoria is disconnected from the rest of 
other vertices of U in HRRNG, HnBS(1.5), and HRGG. Of course, this 
cluster represents the most important economic activity of the region. But in 
this case the connectivity result is poor with respect to intersection with P. A 
similar situation takes place with Merida-Cancun connection in intersection 
of minimum spanning tree MST with highway graph H. 



Let a be a transport link connecting San Luis Potosi region with Mexico city- 
Puebla-State of Mexico region then we see that links of motorway network H 
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which are also edges of Gabriel graph GG correspond to the common links 
of the motorway network H and Physarum graph P(0.19). San Luis Potosi 
connection represents a special case linked directly to Mexico city region with 
motorway M57. 

Finding 7 H R GG [j{a} = H R P(0.19) 

Unfortunately, we have lost a link between Merida-Cancun and Mexico cen- 
ter and also a center-to-north link from San Luis Potosi to Chihuahua. This 
means that motorways M85 and M190 are partly misrepresented. Said that, 
the pacific-long highway connecting Nogales to Mexico city by M15 is pre- 
served. Thus, at least, we have preserved one of the largest and most important 
connection between Mexico city and the North of Mexico. This way, P(0.58) 
is the best approximation as we saw previously. 



5 Discussion 



To approximate, or rather re-construct, development of transport network in 
Mexico we cut of Mexico-shaped plate of agar, represented nineteen major ur- 
ban regions by oat fiakes and placed a Plasmodium of Physarum polycephalum 
in place of Mexico cityp] The Plasmodium developed into a fuU-fiedged Plas- 
modium spanning all, or almost all, oat fiakes on the agar plate. We compared 
graphs of protoplasmic networks developed by Physarum with man-made Fed- 
eral Highway network and with basic types of proximity graphs. 

We found that Physarum-made graph is a sub-graph of the highway network 
apart of few ages. This means that, in principle, slime mould based approx- 
imation of man-made transport networks works even in very simple experi- 
mental setups. The slight mismatch between Physarum and highway graphs 
may be because we did not take geographical profile of Mexico into account 
completely, no mountains or rivers were mapped onto the Plasmodium growth 
substrate. Also, inoculation of Plasmodium in Mexico city was backed up only 
by the fact that Mexico city is most populated region of the country. Thus, 
our future work might focus on experiments with three-dimensional substrates, 
which adequately refiect true landscape of Mexico, and on adding a historical 
perspective to our studies, particularly in choosing more ancient urban region 
for Plasmodium inoculation. 

Comparison of Physarum-made transport networks developed in United King- 
dom [9j to the networks grown in Mexico brings interesting insights on the 



^ You can see an approximation of these simulations from http: //www.youtube . 
fcom/watch?v=0mwtPR2cV-4, 
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existing motorway structure in these countries. Physarum graphs, which built 
of only edges with highest probability of occurrence in experiments were min- 
imum (for UK) and almost-minimum (for Mexico) spanning trees. Physarum 
graph is a sub-graph of UK motorways graph, while Physarum graph is a 
super-graph of Mexican highway graph. Does it mean that Mexican Federal 
Highway are redundant? Or do they simply reflect harsher geographical condi- 
tions which lead to a need of duplicating traffic routes? Further studies might 
bring answers to these questions. 
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